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Abstract 

We introduce four-point functions in the hadronic ladder resummation approach to 
large N c QCD Green functions. We determine the relevant one to calculate the Bk 
kaon parameter in the chiral limit. This four-point function contains both the large 
momenta QCD OPE and the small momenta ChPT at NLO limits, analytically. We 
get B K = 0.38 ± 0.15. We also give the ChPT result at NLO for the relevant four- 
point function to calculate Bk outside the chiral limit, while the leading QCD OPE 
is the same as the chiral limit one. 



1 Introduction 



Indirect kaon CP-violation in the Standard Model (SM) is proportional to the matrix 
element 

(K°\K°) = -iC AS =2C(v)(KV\j ' ^yQ±s=2{y)\K Q ) 

1 fi - 

= -iC± s = 2 -B K Flm\ (1) 

with 

Qas=2(x) = AL»(x)L^x) ■ 2L M (x) = [s lfl (l - 76 )d] (x) . (2) 

and C{v) a Wilson coefficient. C{y) and Cas=2 are known in perturbative QCD to next- 
to- leading order (NLO). A review can be found in PQ. 

The kaon bag parameter, B k, defined in ([I]) is an important input for the unitarity 
triangle analysis and its calculation has been addressed many times in the past. There 
have been four main QCD-based techniques used to calculate it: QCD-Hadronic Duality 
[21 E], three-point function QCD Sum Rules jlj, lattice QCD and the 1/N C (N c = number 
of colors) expansion. Recent reviews of the unitarity triangle, where the relevant references 
for the inputs can be found, are j^j. For recent advances using lattice QCD see jUllZllE]- 

Here, we present a determination of the Bx parameter at NLO in the 1/N C expansion. 
That the l/N c expansion introduced in OHO] would be useful in this regard was first sug- 
gested by Bardeen, Buras and Gerard ^1] and reviewed by Bardeen in [T2"l ITS] . There one 
can find most of the references to previous work and applications of this non-perturbative 
technique. 

Work directly related to this work can be found in flJEI] where a NLO in 1/N C calcu- 
lation of Bk within and outside the chiral limit was presented. There, the relevant spectral 
function is calculated at very low energies using Chiral Perturbation Theory (ChPT), at 
intermediate energies with the extended Nambu-Jona-Lasinio (ENJL) model JH] and at 
very large energies with the operator product expansion (OPE). The method, including 
how to deal with the short- distance scheme dependence, has been discussed extensively 
in The improvement in the present work is basically in the intermediate energies and 
the numerical quality of the matching to short distances. 

Another calculation of Bk in the chiral limit at NLO in the 1/N C expansion is in [T7j . 
There, the relevant spectral function is saturated by the pion pole and the first rho meson 
resonance -minimal hadronic approximation (MHA). In JH], the same technique as in 
[TT] was used but including also the effects of dimension eight operators in the OPE of 
the AS* = 2 Green's function and adding the first scalar meson resonance to the relevant 
spectral function. This works differs from that in two aspects. We use our A-boson method 
which allows for an easy identification of the precise scheme used in the hadronic picture by 
only using currents and densities directly in four dimensions rather than using dimensional 
regularization throughout. This is a difference in formalism, not in physical content. On 
the other hand, we use a much more elaborate way to include hadronic states allowing 
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the use of more constraints JH] and thus more resonances than the work of fT7\ ITS] . The 
hadronic model developed in jTH] also includes current quark masses effects and will be 
used to determine Bk in the real case [20] . 

The present manuscript is organized as follows. In Section El we collect a few perturba- 
tive formulas at NLO in QCD for completeness. Section El explains in detail how we derive 
the value of Bk from a large N c Green's function with two densities and two currents. 
This can be found essentially in [T5] but we add here an extra QCD calculation of the 
short- distance constraint analogous to the equivalent terms in our work on Q 7 and Q$ |21j . 
QCD is still not solved even at leading order in the 1/N C expansion. Therefore, at present, 
any NLO in 1/N C approach to weak matrix elements needs a large N c hadronic model. 
We present in Section |U our hadronic model based on a severe approximation to a ladder 
resummation of QCD as explained in [JJJ]. This model is used at intermediate distances 
and is the main uncertainty in our approach. Then we discuss the chiral limit results in 
Section The extension beyond the chiral limit needs more work both at intermediate 
distances but we include some comments already in Section El At short distances, the 
leading dimension six OPE operator is known and agrees with the chiral one. The last 
section contains our conclusions and a comparison with other approaches. 

Some preliminary results were already presented in and |23j. 



2 Some perturbative QCD formulas 



This section is a compilation of the perturbative QCD coefficients at NLO from pQ needed 
in our analysis of Bk and derived in |24j . 

The coefficients C{y) and Cas=2 are known in perturbative QCD at next-to-leading 
order (NLO) in a = as/ tt in two schemes the 't Hooft-Veltman (HV) scheme (MS 

subtraction and non-anti-commuting 75 in D 7^ 4) and in the Naive Dimensional Regu- 
larization (NDR) scheme (MS subtraction and anti-commuting 75 in D 7^ 4). Cas=2 col- 
lects known functions of the integrated out heavy particle masses and Cabibbo-Kobayashi- 
Maskawa matrix elements, it can be found in ^ . Its actual form is very important for 
CKM analysis but not needed here. 

The Wilson coefficient C{u) is 

72 P2I1 



C{v) 



1 +a(v) 



A P. 

71 is the one-loop AS = 2 anomalous dimension 
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Pi and P2 are the first two coefficients of the QCD beta function 
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The explicit numbers are for N c 



n 



f 



3 X-boson Method and Known Constraints 

The X-boson method was explained in detail in ^H] ■ It takes the idea of ^T] of reducing 
the four-quark operator in (j2J) to products of currents and follows through the full scheme 
and scale dependence. In ^3] it was explicitly showed how short-distance scale and scheme 
dependences can be taken into account analytically in the 1/N C expansion. Here, we only 
sketch the procedure introducing the notation. 
The effective action Ta5=2, 



AS=2 



-C AS= 2 C(u) I d A y Q AS =2(y) + h.c. , (7) 



contains all the short-distance physics of the SM. We replace it by the exchange of a 
colorless heavy AS" = 2 X-boson with couplings 

r LD = 2g AS=2 ( f i c ,---) J d 4 yX^y)L,(y) + h.c. . (8) 

The coupling g A s=2{.Hc, 4 • •) is obtained [TH] with an analytical short-distance matching us- 
ing perturbative QCD. 1 Afterwards we only need to identify the current L M in the hadronic 
picture, not the four-quark operator Qas=2- 
The matching leads to 



9as=2(^c, ■ ■ ■) _ r n( s 
^ = Cas=2C(u) 



1 + a (71 log (—^) + Ar ) 



(9) 



The one-loop finite term Ar is scheme dependent 



Ar NDR = _H f x _ ±\ = _H • Ar HV = _ 7 _( 1 _1_\ (1Q) 

8 V NJ 12 ' 8 V NJ K ' 

and makes the coupling |«7as=2| scheme independent to order a 2 . This coupling is also 
independent of the scale v to the same order. Notice that there is no dependence on the 



1 At energies small compared to the W^-boson mass but where perturbative QCD is still valid. 
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cut-off scale \ic -this feature is general of four-point functions which are product of con- 
served currents [T3] . This procedure thus includes the standard leading and next-to- leading 
resummation to all orders of the large logs in [as log(Mw/v)] n and as [as log(Mjy/V)] n in- 
cluding the short-distance scheme dependence. 

In order to get at the matrix-element (0), we calculate a two-point Green function in the 
presence of the weak effective action, with pseudo-scalar densities carrying kaon quantum 
numbers. After reducing the kaon two-quark densities, the two-point function (fTT|) provides 
the matrix element in (P) via standard LSZ reduction. The two-point function is evaluated 
using the matching with the X-boson effective action. We thus want to calculate the 
two-point function [T^ ITK] 

n A 5= 2 (g 2 ) = ifd i e^(0\T(P^(0)P K a(x)e <r ™) |0> (11) 
which we need to order g\ s =2- After reducing, we obtain 

(X%)|e ir — \K°(q)) = (X%)|e* r ">|X°(g)> = -iC AS =2 y B K q 2 F 2 K 

_ f d 4 p x gjs=2 % 9»v T\^(r? a 2 ) (12) 



where q 2 is the external momentum carried by the kaons. The basic object is the reduced 
four-point function 

n^(p 2 x ,q 2 )=t 2 4(K°(q)\ J d 4 x J d 4 y e ' ipx < x -^ T (L»(x) L v (y)) \K°(q)) . (13) 

This can be obtained from the four-point Green's function with two kaon pseudo-scalar 
densities and two currents L M . 

At large N c , the reduced four-point function (JT5j) factorizes into two disconnected two- 
point functions at all orders in quark masses and external momentum q 2 . This disconnected 
part is 

g^ZconAP^Q 2 ) = (2nrS^(p x )8q 2 F 2 K , (14) 
which leads to the well-known large- iV c prediction 2 

B N k=\. (15) 



At next-to-leading order in the 1/N C expansion, one has 



d 3 g&s=2 



4M|C AS=2 

with Q 2 the X-boson momentum in Euclidean space and 



(16) 



2 In the strict large N c limit we also have C{v) = 1. 
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The next point is the calculation of (|17j). There are two energy regimes where we know 
how to calculate F(Q 2 ) within QCD, namely, at very large Q 2 and at very small Q 2 . 

In the first regime Q 2 >> 1 GeV 2 , with q 2 kept small and we can use the operator 
product expansion in QCD. One gets 



c£ +2 (v,Q 2 )(K"(q)\Q% +2 \K (q)) 



(0 



n=2 i=l 



Q 



2n 



where Q2I+2 are local AS = 2 operators of dimension 2n + 2. In particular, 



Q 6 



(TxL ll (x)L fl (x) 



and 



C 6 (u,Q 2 



^n 2 jia 



1 + a 



(Pi -7i) log( 



wi 



ith 3 



Fi 



72 
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+ (A-7i) 
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7i 



+ Fi 



119 
~16~ 



+ 0(a 2 



(18) 



(19) 



(20) 



(21) 



where the term with Fi was not known before. The finite term Fi is order N c and therefore 
this a 2 term is of the same order in N c as the leading term. In fact, at the same order in 
N c , there is an infinite series in powers of a. 

The above can be used to take the limit Mx 00 explicitly via 



B 



K 



-C(u) ( 1 + a ( 7l log ( ^ ) + Ar 



v 



16n 2 F 2 



dQ 2 F[Q 



2]M X -^00 



+ 



dQ 2 F[Q 2 )^ 8 



2]M x ^oo 



+ 



Ml, 



+ 0{a 2 



(22) 



Where -F[Q 2 ]d>8 is obtained inserting in (JT7j) the result in ()18p minus the dimension six 
term. 

For the list and a discussion of the dimension eight operators see |18[l2f)|. In ^Hj there is 
a calculation in the factorizable limit of the contribution of the dimension eight operators. 
Numerically, the finite term of order a 2 competes with that contribution when Q 2 is around 
(1 ~ 2) GeV 2 . 

The second energy regime where we can calculate -F[Q 2 ] model independently is for 
Q 2 — ► 0, where the effective quantum field theory of QCD is chiral perturbation theory. In 
ChPT the result is known up to order p 4 both in the chiral limit ^3 Ejj and outside the 
chiral limit |22|. The result in the chiral limit is 



F*[Q 2 ] = 3 + A 4 Q 2 + A 6 Q 4 



12 

A 4 = --2 (2L X + 5L 2 + L 3 

-Tn 



Lc 



(23) 



3 The leading terms differs from the one in ^7] only in terms subleading in 1/N C 
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with Fq the chiral limit of the pion decay constant F n = 92.4 MeV. The next term, Aq, 
can be easily calculated but contains some of the unknown constants from the p 6 ChPT 
Lagrangian. 

We still need to describe the intermediate energy region for which we use the large N c 
hadronic model described in the next section. We have used that model to predict the 
series in ()23j) minus 3 + A^Q 2 , which is known. This is equivalent to predict the relevant 
Ci and higher couplings combinations. 

4 A Ladder Resummation Large N c Model 

All the results and information we have presented so far on B k are model independent. In 
particular, we have seen in the previous section that there are two energy regimes which 
can be calculated within QCD. In this section, we describe a large N c hadronic model that 
provides the full n^ nn {Q 2 ,q 2 ) which, apart from other QCD information, contains these 
two QCD regimes analytically. 

The large N c hadronic model we use was introduced in [T^j. It can be thought of as 
QCD in the rainbow or ladder-resummation approximation. The basic objects are vertex 
functions with one, two, three, . . . two-quark currents or density sources attached to them, 
referred to as one-point, two-point, three-point, . . . vertex functions. These correspond to 
the two-particle irreducible diagrams in large N c QCD. These vertex functions are glued 




Figure 1: One-point, two-point, three-point, • ■ ■ vertex functions. The crosses can be 
vector or axial-vector currents and scalar or pseudo-scalar densities. 

into infinite geometrical series with couplings gy for vector or axial-vector sources and g$ 
for scalar or pseudo-scalar sources, what can be seen as a very crude approximation to the 
two-particle reducible part. In this way one can construct full n-point Green's functions in 
the presence of current quark masses -see for instance, how to get full two-point functions 
in Figure El 




Figure 2: Infinite geometrical series which gives full two-point functions at large N c . The 
black vertices that glue the vertex functions together are either gy for vector or axial-vector 
sources or g$ for scalar and pseudo-scalar sources. 
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The basic vertex functions in Figure ^ have to be polynomials in momenta and quark 
masses to keep the large N c structure. As explained in JH] , only the first nonet of hadronic 
states per channel, i.e., the pseudo-scalar pseudo-Goldstone bosons, the first vectors, the 
first axial-vectors and the first scalars can be easily generated within this framework. The 
coefficients of the vertex functions are free constants of order N c , many of which can be 
fixed by imposing chiral Ward identities on the full Green's functions. Chiral perturbation 
theory at order p 4 and the operator product expansion in QCD help to determine many 
more of these free coefficients in the vertex functions. A consequence of the formalism 
described here is that the vertex functions obey Ward identities with a constituent quark 
mass [TH] . 

The full two-point Green's functions obtained from the resummation in Figure El agree 
with the ones of large N c QCD when one limits the hadronic content to be just one hadronic 
state per channel -our model does not produce any less or more constraints than large N c 
QCD and all parameters can be fixed in terms of resonance masses i n agreement with 
other groups using large N c and one state per channel. Introducing two or more hadronic 
states per channel systematically is difficult as explained in We leave for future work 
the investigation into how to carry it out. 

Most low-energy hadronic effective actions used for large N c phenomenology are in the 
approximation of keeping the resonances below some hadronic scale and always in the 
chiral limit. In many cases less than the four states included here are used. 

The procedure described above of obtaining full Green's functions can be done in the 
presence of current quark masses. In fact, in (T^j two-point functions were calculated 
outside the chiral limit and all the new parameters that appear up to order m 2 q can be 
determined except one; namely, the second derivative of the quark condensate with respect 
to quark masses. Some predictions of the model we are discussing involving coupling 
constants and masses of vectors and axial-vectors in the presence of masses are 



where i,j, k, I are indices for the up, down and strange quark flavors. 

Three-point functions, as shown in Figure El were calculated in the chiral limit in ^H] ■ 
We have now all the needed ones in the study of Bk also outside the chiral limit and they 
will be presented elsewhere j2E]. Some three-point functions have also been calculated in 
other large N c approaches and in the chiral limit j23 EU I2H EDI ~for instance, PVV, PVA, 
PPV and PSP three-point functions, where P stands for pseudo-scalar, V for vector, A for 
axial-vector and S for scalar sources. They agree fully with the ones we get in our model 
when restricted to just one hadronic state per channel. 

Four-point functions are constructed analogously with the two-topologies dictated by 
large N c , one with two three-point vertex functions and one with a four-point vertex 




fvki M Vk i , 

flkl ^Akl + fkl > 

~(qq) x (raj + m,- -m k - mi) , 



(24) 
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V 

Figure 3: Infinite geometrical series which gives full three-point functions at large N c . The 
crosses that glue the vertex functions are either gy for vector or axial-vector sources or g$ 
for scalar and pseudo-scalar sources. 

function. The final four-point functions fulfill the correct chiral Ward identities -including 
current quark masses and we impose short-distance QCD constraints as well as the ChPT 
results at NLO to determine the free parameters allowed by the chiral Ward identities. 



5 Chiral Limit Results 

We now calculate -F[Q 2 ] of Eq. (|17|) in the chiral limit using the model of the previous 
section. After integrating over the four-dimensional Euclidean solid angle Qq and doing 
the limit q 2 — > as in ([T7]). we get 



; -\ vVi , /Y> , a V OCA as Pv 



+ ™ + JX + 'A (25) 



13a IV lA 

(Q* + M\f (Q^ + m v f (Q* + Mi) 

We have used here the two-point vertex functions to second order in momenta, higher 
orders will make it impossible to satisfy the Weinberg sum rules. The three- and four- 
point vertex functions have been expanded to fourth order in the external momenta. The 
short- distance constraints used are the Weinberg sum rules and the fact that the vector 
form factor should vanish as 1/Q 2 for large Q 2 . As expected from the arguments in [T5] . 
we find clashes between short- distance constraints for three-point functions and those of 
the needed four-point function here. I.e., with just a finite number of hadronic states per 
channel not all short-distance constraints for n-point Green are compatible in general. Of 
course, one can always impose the latter which are compatible with the subset of the short- 
distance constraints corresponding to the momenta structure of the three-point functions 
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— 9 91 


10 3 L 5 


0.97 


1.46 


10 3 L 9 


5.93 


6.9 


10 3 L 10 


-4.4 


-5.5 


My 


805 MeV 


690 MeV 


M A 


1.16 GeV 


895 MeV 


M s 


1.41 GeV 


1.06 GeV 


A, 


-12.7 GeV~ 2 


-23 GeV" 2 



Table 1: The low energy inputs used together with some derived quantities for the two sets 
of input. 

entering in the four-point functions but not with all three-point functions short- distance 
constraints. That is what we have done for calculating Bk- 

The short-distance constraints on -F[Q 2 ] discussed in Sect. El and Ref. |TZl HHj require 
Si = 5 2 = in (j23j). Imposing these we obtain the expression for F[Q 2 ] given in the 
appendix. 

The explicit calculation reveals that diagrams with the exchange of vector and axial- 
vector states produce not only single poles but also double and triple poles. These also 
arise from diagrams with only one of these propagators due to the factor 1/q 2 present in 

(Hlj). 4 

The function F X [Q 2 } in (|T7|) reproduces the large iV c -pole structure found in jT?l ITHJ 
but including the first hadronic state in all the spin zero and one channels. 

The relevant free parameters are the masses My, Ma and Ms, the pion decay constant 
in the chiral limit, F Q , and three combinations of constants appearing in the three- and 
four-point vertex functions, ^1,2,3- 

The low energy inputs we use are Fq and the 0(p 4 ) couplings Li determined from the 
ChPT fits to data. We use two different fits, namely, fit 10 of [3T] for £1,2,3 and L910 fit 
from both with the full p 6 fits and the ones using only p 4 expressions. The input 
values for these two cases p 6 and p 4 are given together with some derived quantities in 
Table [T] The values of the masses and the slope follow from the Li and F used as 
input using the relations of [TH] and (J2"3~j) . One could also use the physical masses and then 
use the Li only to get at the slope. The case with p e input is within 30% of the physical 
masses. For the scalars, which mass to use is still an open question, but with the mounting 
evidence that the first one are not present in large N c , and references therein, a mass 
of around 1.4 GeV seems fine. 

4 Ref. |171 118| use a different underlying four-point function. The same pole structure does show up 
there as well as is required by chiral symmetry. 
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Q 2 [GeV 2 ] Q 2 [GeV 2 ] 

(a) (b) 

Figure 4: F[Q 2 ] from ChPT, short distance and the large N c ladder resummation model 
interpolation, (a) with the p 6 fit values as input (b) with the p 4 fit values as input. 

We determine A\ from the slope v4 4 as given in (J23*j) and the numerical values from 
Tabled 

A% and A% are obtained using the short distance constraints of ()18j) . via 

F[Q 2 ]sd = p + £ + O (£) . (26) 

We use the estimate D 8 /D e ps 0.13 GeV -2 of [IH] but with rather large uncertainties. 
is calculated from (J2*U|) which contains the value of i?x inside it. One can choose here the 
large N c values or the values which come out of the analysis. We choose the latter ones 
using Bjc = 0.17,0.37 for the p 6 ,p 4 inputs, respectively. With a s (m T ) = 0.35 we obtain 
D 6 = 0.0.028, 0.052 GeV 2 and D 8 = 0.0036, 0.0068 GeV 4 . In the remainder we use the 
short- distance (|2l)j) with these values. 

We have then used two different ways to match the model F[Q 2 ] to short distance. The 
first one is to have the two coincide at the values of Q 2 = 2 and 3 GeV 2 , these are labeled 
with "values" in the figures. The other is to have the model reproduce the values of D e 
and D 8 in its large Q 2 expansion, this case is labeled "order" in the figures. 

In Fig. 0] we have plotted the short-distance curve and the model curves as well as the 
order p 4 ChPT approximation. 

We now use (J22|) and integrate (fTTj) up to the matching point // and from that point 
on we use the OPE result including dimension eight corrections ^H] of (|26|. This OPE 
contribution to B\ is negligible. In Figure we plot this result as a function of fi 2 Notice 
the nice plateau one gets between 1 and 2 GeV 2 , with the exception of the case "order" 
with the p G input. From FigHk it can be seen that this because this case has an extremely 
slow approach to the short- distance. Varying the inputs within does not change the result 
more than those shown. We thus obtain 

B\ = 0.38 ±0.15 (27) 
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Figure 5: B\ plotted vs the upper limit of the integral in p7jl . See text for further 
explanation. 

which is fully compatible with the one found previously also in the chiral limit in [T4*l \W\ 

Ed EE]. 

6 Outside the Chiral Limit Results 



The QCD results of the two energy regimes explained in Section |3] are known outside the 
chiral limit too. 

The function -F[Q 2 ] is known for very large values of Q 2 , which can be calculated using 
the OPE in QCD. In fact, the dimension six operator and Wilson coefficient are the same 
as the chiral limit ones. Differences will first start to appear at dimension 8 and we expect 
them to be small given the overall small estimate of the dimension 8 corrections. 

At small values of Q 2 we can use ChPT. It is easy to get that F[0] = for the real 
case instead of F x [0] = 3 in the chiral limit. Notice that the value of F x [0] is a strong 
constraint on the value of B\. That F[0] = outside the chiral limit is a strong indication 
that corrections to the leading in N c value E>k = 3/4 will be much smaller than in the 
chiral limit case. 

The ChPT calculation in the real quark masses case of -F[Q 2 ] has been done to order 
p 4 [22] and -F[Q 2 ] remains small -below 0.15- up to energies around 0.2 GeV 2 and then 
goes negative. 



Outside the chiral limit and to chiral order p , we get 



21(2) 



Q 6 + 2m\Q 4 + 2m 4 K Q 2 



m 



K 



+ m 2 K ) 2 



+ 4 



1 Q 2 + 2ml + 2m 2 



+ 4 



3 Q 2 + 2m 2 K + 2m 2 



in 



i[Q 



2 2 

m K , 



in 



K 



[2%) 
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To chiral order p A and leading in 1/N C we get 



F[Q 2 } 




(29) 



The angular integration results in the function 



I[x,m 1 ,m 2 ] 



x — mi + ^2 



1 + Ami x/(x — mi + — 1 



(30) 



Inii 



We have used the Gell-Mann-Okubo relation to simplify the p 4 contribution and dropped 
all terms that are subleading in 1/N C as well as changed _F 2 in the higher order to F 2 . 

In Fig. El we have shown -F[Q 2 ] at low Q 2 from the formulas above in ChPT for defi- 
niteness with the p 6 inputs of Table H It can be seen that the function starts at zero and 
the area under the curves, which gives the corrections is much smaller away from the chiral 
limit. The discontinuity in the curve is due to the fact that to the right the state with an 
intermediate pion can go on-shell. 

One expects that higher ChPT order terms correct the behaviour of -F[Q 2 ] for Q 2 > 
0.2 GeV 2 and -F[Q 2 ] will tend to the chiral limit curve in Figure 01 The question is at 
which energy does it happens. This can only be answered using a hadronic large N c 
approximation to QCD at present. So, though we have strong indications that the value 
Bk = 3/4 has small chiral corrections as shown before, we have to wait till we get the full 
four-point Green's function in the real case (|13|) to confirm it. The number of new free 
parameters in our hadronic model including quark masses is rather high and we have at 
present not been able to determine all as we did in the chiral limit discussed earlier. We 
will eventually present the full result in the presence of current quark masses in [2*Uj . 

7 Summary and Conclusions 

In the past few years, lattice QCD has produced many calculations of Bk using different 
fermion formulations -see jEl 13 E] for references. Those include some chiral limit extrap- 
olations -like for instance in [H3] obtaining B K = 0.32 ± 0.22 with quenched staggered 
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Figure 6: The ChPT calculation of -F[Q 2 ] in and away from the chiral limit. 

quarks, in obtaining B\ = 0.34 ±0.02 with two dynamical domain wall quark flavours 
or in [36J obtaining B\ = 0.39 ± 0.03 with quenched domain wall quarks -which show a 
clear decreasing tendency with respect to the real case value, in agreement with the results 
found here and in |TH El El E|- Promising preliminary unquenched results have also 
started to appear in the last two years [331 EH] • 

The QCD-Hadronic Duality result for -B^|21E1 is very close to the chiral limit result 
above because -as already mentioned in what was calculated there, is the order p 2 
coefficient of the chiral expansion which actually is the chiral limit value of Bk- 

We have found values in this paper which are very compatible with the those chiral limit 
values and in agreement with those of El El E] ■ The improvement over our earlier 
work [T^ ITo] is the use of a more reliable model for the intermediate energy regime. What 
we have done beyond the work of ^3 E] is to include more resonances and a somewhat 
different input for the low-energy constants. In addition we have argued that away from 
the chiral limit the corrections to the LO large N c value 3/4 should be much smaller. A 
full result in the presence of current quark masses will be presented in |20j . 
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A The Explicit Form of F[Q 2 } 

The function -F[Q 2 ] has in the chiral limit the form of Eq. (|25J) . The short distance already 
imposes 

5 1 = 5 2 = 0. (31) 

The remaining constants are 

QM V 3M v 3 2 15 My 12M V 3 6 

aV ~ ~ (M v - M-f ~ My— Ml ~ 4 V + TM* + -MT A * + 128 M ^ 

6M V 3M V 3M v Mg ^_^A_^ M 2 9M v 9 M 2 

aA ~ (My - Ml) 2 + M v -M\~ ( M \ - Ml) 2 + " T My + AM\ ~ 2 

12M| 12MJ „ _ e 

Ml 1 (Ml- Ml) 2 2 3 2 V 3 384M 2 



3M v Ml + Y1M%A\ 
(Ml - Ml) 2 



3M v Ml 9 3M V 6Mf 3 

PV ~ ~(M V - Ml) 2 + l Mv ~^l' ~Ml A " ~ 64 y 



3M V 9M V 3M v Ml 15 Ml 3 „ r2 „ r2 3 „ r , . , r 



y 1 Ml -Ml 2 v a s 192 My 
3M v Ml 1 MA 2 

7A = - Mv -Ml + 3MMf N « (32) 
The parts proportional to N n defined by 

N-^-(l-W\ (33) 

come from the diagrams with anomalous three-point vertex functions. 
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